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Abstract 

Wc numerically construct asymptotically Anti-dc Sitter charged black hole so- 
lutions of (4 + l)-dimensional SU (2) Einstein- Yang-Mills theory that, for suffi- 
ciently low temperature, develop vector hair. Via gauge-gravity duality, these 
solutions describe a strongly-coupled conformal field theory at finite tempera- 
ture and density that undergoes a phase transition to a superfluid state with 
spontaneously broken rotational symmetry (a p-wave superfluid state). The 
bulk theory has a single free parameter, the ratio of the five-dimensional grav- 
itational constant to the Yang-Mills coupling, which we denote as a. Previous 
analyses have shown that in the so-called probe limit, where a goes to zero 
and hence the gauge fields are ignored in Einstein's equation, the transition to 
the superfluid state is second order. We construct fully back-reacted solutions, 
where a is flnite and the gauge fields are included in Einstein's equation, and 
find that for values of a above a critical value = 0.365 ± 0.001 in units of the 
AdS radius, the transition becomes first order. 
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1. Introduction 

The Anti-de Sitter/Conformal Field Theory correspondence (AdS/CFT) [1] 
provides a novel method for studying strongly-coupled systems at flnite density. 
As such, it may have useful applications in condensed matter physics, espe- 
cially for studying scale-invariant strongly-coupled systems, for example, low- 
temperature systems near quantum criticality (see for example refs.[2, 3, 4] and 
references therein). Such systems are not purely theoretical: the (thcrmo)dynamics 
of some high-Tc superconductors may be controlled by a quantum critical point. 
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AdS/CFT is a holographic duahty: it equates a weakly-coupled theory of 
gravity in {d+ l)-dimensional AdS space with a strongly-coupled d-dimensional 
CFT "living" at the AdS boundary. CFT states with finite temperature are 
dual to black hole geometries, where the Hawking temperature of the black 
hole is identified with the temperature in the CFT [5] . The current of a global 
U{1) symmetry in the CFT will be dual to a U{1) gauge field in AdS space. 
AdS/CFT thus allows us to compute observablcs in a strongly-coupled CFT, 
in states with finite temperature and density, by studying asymptotically AdS 
charged black holes. By now AdS/CFT can model many basic phenomena in 
condensed matter physics, such as the quantum Hall effect [6], non-relativistic 
scale-invariance [7, 8], and Fermi surfaces [9, 10, 11]. 

AdS/CFT can also describe phase transitions to superfluid states, that is, 
phase transitions in which a sufficiently large U{1) charge density triggers spon- 
taneous breaking of the U{1) symmetry: an operator charged under the U{1) 
acquires a nonzero expectation value [12, 13, 14]. We will refer to this as the 
operator "condensing." The simplest bulk action that can describe such a tran- 
sition is Einstein-Maxwell theory coupled to a charged scalar. In the bulk, a 
charged black hole develops scalar hair. In the CFT, a charged scalar operator 
condenses. 

A simple bulk action has one great virtue, namely a kind of universality: the 
results may be true for many different dual CFT's, independent of the details 
of their dynamics. For the Einstein-MaxwcU-scalar case, a fruitful exercise is to 
study various functional forms for the scalar potential and to scan through values 
of couplings in that potential [15, 16, 17]. Generally speaking, scanning through 
values of these parameters corresponds to scanning through many different dual 
CFT's. As shown in refs. [15, 16, 17], such changes can have a dramatic effect, 
for example, the phase transition can change from second to first order. 

AdS/CFT can also describe superfluid states in which the condensing op- 
erator is a vector and hence rotational symmetry is broken, that is, p-wave 
superfluid states [18, 19]. Here the CFT has a global SU{2) symmetry and 
hence three conserved currents J^, where a = 1,2,3 label the generators of 
SU{2). For a sufficiently large charge density for some U{1) subgroup of SU{2), 
say a sufficiently large (J3), holographic calculations reveal that, of the known 
available states, those with lowest free energy have a nonzero (Jf). Not only 
is the U{1) broken, but spatial rotational symmetry is also broken to some 
subgroup. 

On the AdS side, a simple bulk action that can describe such a transition is 
Einstein- Yang-Mills theory with gauge group SU{2). CFT states with nonzero 
(J3) are dual to black hole solutions with nonzero vector field Af{r) in the bulk, 
where r is the radial coordinate of AdS space. States with nonzero ( Jf ) are dual 
to black hole solutions with a nontrivial A], (r) . The superfluid phase transition 
is thus dual to charged AdS black holes developing vector hair. A string theory 
realization for this model is given in refs. [20, 21, 22, 23]. 

Unlike the Einstein-MaxwcU-scalar case, SU{2) Einstein- Yang-Mills theory 
has only a single free parameter, a = AC5/5, where K5 is the gravitational con- 
stant (we will work in (4 + 1) dimensions, hence the subscript) and g is the 
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Yang-Mills coupling. The Yang-Mills source terms on the right-hand-side of 
Einstein's equation are proportional to c? . To date, most analyses of the holo- 
graphic p-wave superfluid transition have employed the so-called probe limit, 
which consists of taking a so that the gauge fields have no effect on the 
geometry, which becomes simply AdS-Schwarzschild. The probe limit was suf- 
ficient to show that a superfluid phase transition occurs and is second order. 

Our goal is to study the effect of finite a, that is, to study the back-reaction 
of the gauge fields on the metric. We will work with (4 -|- l)-dimensional SU (2) 
Einstein- Yang-Mills theory, with finite a. We will numerically construct asymp- 
totically AdS charged black hole solutions with vector hair (for similar studies 
see refs. [24, 25, 26]). Our principal result is that for a sufficiently large value of 
a the phase transition becomes first order. More specifically, we find a critical 
value ttc = 0.365 ± 0.001 in units of the AdS radius, such that the transition is 
second order when a < ac and first order when a > ac- 

We can provide some intuition for what increasing a means, in CFT terms, 
as follows. Generically, in AdS/CFT I/k^ oc c, where c is the central charge of 
the CFT [27, 28], which, roughly speaking, counts the total number of degrees 
of freedom in the CFT. Correlation functions involving the SU{2) current will 
generically be proportional to [2. 3, 4]. We may (again roughly) think of 
l/g"^ as counting the number of degrees of freedom in the CFT that carry SU (2) 
charge. Intuitively, then, in the CFT increasing a means increasing the ratio of 
charged degrees of freedom to total degrees of freedom. 

The paper is organized as follows. In section 2 we write the action of our 
model and discuss our ansatz for the bulk fields. In section 3 we describe 
how to extract thermodynamic information from our solutions. In section 4 we 
present numerical results demonstrating that increasing a changes the order of 
the phase transition. We conclude with some discussion and suggestions for 
future research in section 5. 



2. Holographic Setup 

We consider SU{2) Einstein- Yang-Mills theory in (4-|-l)-dimensional asymp- 
totically AdS space. The action is 



hdy • 



(1) 



where K5 is the five-dimensional gravitational constant, A = — is the cosmo- 
logical constant, with L being the AdS radius, and g is the Yang- Mills coupling 
constant. The SU{2) field strength F^^ is 



= d^Al - d^Al + e'^^'^A^Al , (2) 

where /U, = {t,r,x,y, z}, with r being the AdS radial coordinate, and e"*"^ is 
the totally antisymmetric tensor with e^"^^ = +1. The A'Jj^ are the components 
of the matrix-valued gauge field, A = A'^r'^dx'^, where the r" are the 5/7(2) 
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generators, which are related to the Pauli matrices by t° = a"' /2i. Sbdy inchidcs 
boundary terms that do not affect the equations of motion, namely the Gibbons- 
Hawking boundary term as well as counterterms required for the on-shell action 
to be finite. We will write Sj^dy explicitly in section 3. 

The Einstein and Yang-Mills equations derived from the above action are 

R/xv + j2^i^'^ ~ '^5 ^-^M'' ~ '^Tp^'g^j^J^ , (3) 

where the Yang-Mills stress-energy tensor T^j,^ is 

V = [Pp^F", - -^g^^Fp^F'"'^ . (5) 

Following ref. [18], to construct charged black hole solutions with vector 
hair we choose a gauge field ansatz 

A = 4>{r)T^dt + w{r)T^dx . (6) 

The motivation for this ansatz is as follows. In the field theory we will introduce 
a chemical potential for the U{1) symmetry generated by r^. We will denote 
this U{1) as U{1)3. The bulk operator dual to the C/(l)3 density is A^, hence 
we include A^{r) = (j)(r) in our ansatz. We want to allow for states with a 
nonzero (Jf), so in addition we introduce AjJ,(r) = w{r). With this ansatz 
for the gauge field, the Yang-Mills stress-energy tensor in eq. (5) is diagonal. 
Solutions with nonzero w{r) will preserve only an SO{2) subgroup of the SO{2>) 
rotational symmetry, so our metric ansatz will respect only SO{2). Furthermore, 
given that the Yang-Mills stress-energy tensor is diagonal, a diagonal metric is 
consistent. We will also pattern our metric ansatz after the ones used in ref. 
[26] since these tame singular points in the equations of motion. Our metric 
ansatz is 

ds2 = -N{r)(j{rfAe + ;^dr2 + r'' f{r)-^dx'' + r^f{rf {dy^ + dz^) , (7) 

with N{r) = -^ + ^. For our black hole solutions we will denote the 
position of the horizon as r^. The AdS boundary will be at r — oo. 

Inserting our ansatz into the Einstein and Yang-Mills equations yields five 
equations of motion for m(r), tT(r), /(r), (t>{'r), w{r) and one constraint equation 
from the rr component of the Einstein equations. The dynamical equations can 
be recast as (prime denotes ^) 
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m! 



f" 




6iVcr2 

_a^fw^ ^a'fw'^ _ ^, (i f ^N' (8) 
3r2Ar2o-2 



The equations of motion are invariant under four scaling transformations 
(invariant quantities are not shown), 

(7) a ^ Xa, (j) ^ X(f>, 



{II) f^Xf, w^X 



-2 



W. 



(Ill) r ^ Xr , m ~> A^m , w — > Aw , (p ^ ^4'} 

{IV) r ^ Xr , TO — > A^TO , L — > XL , — > ^ , a — > Xa, 

where in each case A is some real positive number. Using (I) and (II) we can 
set the boundary values of both cr(r) and /(r) to one, so that the metric will 
be asymptotically AdS. We are free to use (III) to set to be one, but we will 
retain rh as a bookkeeping device. We will use (IV) to set the AdS radius L to 
one. 

A known analytic solution of the equations of motion is an asymptotically 
AdS Reissner-Nordstrom black hole, which has (j){r) = fi — q/r^, w{r) = 0, 

a(r) = /(r) = 1, and N{r) ^ (r^ ^ + ^-^), where toq = ^ + and 

q = firl^. Here n is the value of (j){r) at the boundary, which in CFT terms is 

the U{1)3 chemical potential. 

To find solutions with nonzero w{r) we resort to numerics. We will solve 
the equations of motion using a shooting method. We will vary the values of 
functions at the horizon until we find solutions with suitable values at the AdS 
boundary. We thus need the asymptotic form of solutions both near the horizon 
r = rh and near the boundary r = oo. 

Near the horizon, wc define eh = — — 1^1 and then expand every function 
in powers of Ch with some constant coefficients. Two of these we can fix as 
follows. We determine rh by the condition N{rh) = 0, which gives that m{rh) = 
Additionally, we must impose (r^J = (f){rh) = for A to be well-defined 
as a one-form (see for example ref. [29]). The equations of motion then impose 
relations among all the coefficients. A straightforward exercise shows that only 
four coefficients are independent, 

{cl>1,atf^,w^} , (9) 
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where the subscript denotes the order of eh (so ctq is the value of (T(r) at the 
horizon, etc.). All other near- horizon coefficients are determined in terms of 
these four. 

Near the boundary r = oo we define = (^) <^ 1 and then expand every 
function in powers of ei, with some constant coefficients. The equations of motion 
again impose relations among the coefficients. The independent coefiicients are 

{ml^i,4>lwlf^} , (10) 

where here the subscript denotes the power of e;,. All other near-boundary 
coefficients are determined in terms of these. 

We used scaling symnicitrics to set (Jq = /q = 1. Our solutions will also have 
Wq = since we do not want to source the operator Jf in the CFT (J7(l)3 will 
be spontaneously broken). In our shooting method we choose a value of /x and 
then vary the four independent near-horizon coefficients until we find a solution 
which produces the desired value of fj. and has CTq = /o = 1 and = 0. 

In what follows we will often work with dimensionless coefficients by scaling 
out factors of r/j. We thus define the dimensionless fimctions m(r) = m(r)/r^, 
(f){r) = (j){r)/rh and 'w{r) = w{r)/rh, while /(r) and cr(r) are already dimen- 
sionless. 



3. Thermodynamics 

Next we will describe how to extract thermodynamic information from our 
solutions. Our solutions describe thermal equilibrium states in the dual CFT. 
We will work in the grand canonical ensemble, with fixed chemical potential (x. 

We can obtain the temperature and entropy from horizon data. The tem- 
perature T is given by the Hawking temperature of the black hole. 



(11) 




Here k = is the surface gravity of the black hole, with ^ being the 

norm of the timelike Killing vector, and in the second equality we write T in 
terms of near-horizon coefficients. In what follows we will often convert from 
Th to T simply by inverting the above equation. The entropy S is given by the 
Bekenstein-Hawking entropy of the black hole, 

2^ _27ry3_27r^ 3 I2^af 

S - -2^h - —2-rh - -^VT ^ 3 , (12) 

''s ''s '^5 (l2af - (</)^) a2j 

where Ah denotes the area of the horizon and V = J d^x. 

The central quantity in the grand canonical ensemble is the grand potential 
fl. In AdS/CFT we identify fl with T times the on-shell bulk action in Euclidean 
signature. We thus analytically continue to Euclidean signature and compactify 
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the time direction with period l/T. We denote the Euchdean bulk action as / 
and /on-sheii as its on-shell value (and similarly for other on-shell quantities) . Our 
solutions will always be static, hence /on-sheii will always include an integration 
over the time direction, producing a factor of l/T. To simplify expressions, we 
will define / = I/T. Starting now, we will refer to / as the action. / includes a 
bulk term, a Gibbons-Hawking boundary term, and counterterms, 

i = /bulk + /gh + icT ■ (13) 

-^bui^°" and iQ^ff^^"^"^ exhibit divergences, which are canceled by the counterterms 
in /cT- To regulate these divergences we introduce a hypersurface r = r^dy with 
some large but finite rbdy. We will always ultimately remove the regulator by 
taking rb^y oo. Using the equations of motion, for our ansatz /buik^''^' 



iZt^ = ;j«riVa(rV')' • (14) 



=rbdy 

For our ansatz, the Euclidean Gibbons-Hawking term is 



where 7 is the induced metric on the r = rbdy hypersurface and n^dx^ = 
N(r) dr is the outward-pointing normal vector. The only divergence in 
-^bSk'^'*" + /qh'''^''" comes from the infinite volume of the asymptotically AdS 
space, hence, for our ansatz, the only nontrivial counterterm is 

on-shell _ ^ /"iS r- 3 



^CT = / d'^Vr = -^rWNa . (16) 



'5 



=rbdy 



Finally, O is related to the on-shell action, /on-sheih as 

O = lim ion-shell- (17) 

The chemical potential ji is simply the boundary value of Af{r) = (f>{r). The 
charge density ( J|) of the dual field theory can be extracted from /on-sheii by 

(J*) = l lim SF^^-^T^ ,4>l (18) 

Similarly, the current density (Jf ) is 

(Jf) = ^ hm /f°7"V +'^^^ — 3-?- (19) 

V r...^o. 5 Alindy) nl ^^2af_^^Hf^2y' 
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The expectation value of the stress-energy tensor of the CFT is [28, 30] 



Irn \ r ^ ^-^on-shell n. 

(Tij) = hm — ^ .. = hm 



r=rhdy 

(20) 

where i,j,l = {t,x,y,z} and Kij = ^y^WJr) drjij is the extrinsic curvature. 
We find 

''s (^l2af - ((A?*) a2j 
(I^x.) = "^VT^- ,4 K - 8/1) , (21) 



Notice that {Tfx) = {Tty) = {Ttz) = 0. Even in phases where the current 
(Jf ) is nonzero, the fluid will have zero net momentum. Indeed, this result is 
guaranteed by our ansatz for the gauge field which implies a diagonal Yang-Mills 
stress-energy tensor and a diagonal metric (the spacetime is static). 

For m[; = i + <t(1 = 1, = 2/i, /| = 0, and 4)\ = -fi wc recover 

the correct thermodynamic properties of the Reissner-Nordstrom black hole, 
which preserves the 50(3) rotational symmetry. For example, we find that 
(Txx) = {Tyy) = {Tzz) and Vl = —{Tyy). For solutions with nonzero (Jf), the 
50(3) is broken to 50(2). In these cases, we find that {T^x) 7^ {Tyy) = (T^^). 
Just using the equations above, we also find O = —{Tyy). In the superfluid 
phase, both the nonzero ( Jf ) and the stress-energy tensor indicate breaking of 
50(3). 

Tracelessness of the stress-energy tensor (in Lorentzian signature) implies 

{Ttt) = {Txx) + {Tyy) + {Tzz), whlch is indeed true for cq. (21), so in the dual 
CFT we always have a conformal fluid. The only physical parameter in the 
CFT is thus the ratio /x/T. 



4. Phase Transitions 

In this section we present our numerical results. We scanned through values 
of a from a = 0.032 to a = 0.548. Typical solutions for the metric and gauge 
field functions appear in Figure 1. The solutions for other values of a are qual- 
itatively similar. Notice that all boundary conditions are met: at the horizon 
(f){r) vanishes, and at the boundary/Q = ctq = 1 and Wq = 0. 

For every value of a that we use, we find Reissner-Nordstrom solutions for 
all temperatiireS; and for sufficiently low temperatures we always find addi- 
tional solutions, with nonzero w{r), that are thermodynamically preferred to 
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(a) 



(b) 



Figure 1: (a) The dimensionless gauge field components <t>{r) (red dashed) and w{r) (green 
dot-dashed) and the dimensionless metric function rh(r), scaled down by a factor of 10, (black 
solid) versus the AdS radial coordinate r for a = 0.316 at T 0.45rc. (b) The dimension- 
less metric functions o-(r) (red dashed) and /(r) (green dot-dashed) versus the AdS radial 
coordinate r for o = 0.316 at T 0.45Tc. 

the Reissner-Nordstrom solution. In other words, for every value of a that we 
use, we find a phase transition, at some temperature Tc, in which a charged 
black hole grows vector hair, which in the CFT is a p-wave superfluid phase 
transition. Our numerical results show that the phase transition is second order 
for a < ttc and first order for a > ac where Uc ~ 0.365 ± 0.001. 

For example, for a = 0.316 < ac, we only find solutions with (Jf ) = until 
a temperature Tc where a second set of solutions, with nonzero (Jf), appears. 
Figure 2 shows that (Jf ) rises continuously from zero as we decrease T below 
Tc- Figure 3 (a) shows the grand potential V,, divided by n'^VT^/nl, versus 
the rescaled temperature T/Tc for a = 0.316. The blue sohd curve in Figure 3 
(a) comes from solutions with ( Jf ) — and the red dashed curve comes from 
solutions with ( Jf ) ^ 0. We see clearly that at T < Tc the states with ( Jf ) ^ 
have the lower (^tt'^VT^) and hence are thermodynamically preferred. We 

thus conclude that a phase transition occurs at T = T^. The nonzero (Jf) 
indicates spontaneous breaking of ?J(1)3 and of SO(3) rotational symmetry 
down to 50(2), and hence is an order parameter for the transition. Figure 3 (b) 
shows the entropy S, divided by 2tt'^VT^ / k^, versus the rescaled temperature 
T/Tc for a = 0.316. The blue solid curve and the red dashed curve have the same 
meaning as in Figure 3 (a). Here we see that k^S/ (2Tr^VT^) is continuous but 
has a kink, i.e. a discontinuous first derivative, clearly indicating a second-order 
transition. For other values of a < ac, the figures are qualitatively similar. 

A good question concerning these second-order transitions is: what are the 
critical exponents? In the probe limit, a = 0, an analytic solution for the gauge 
fields exists for T near Tc [21], which was used in ref. [19] to show that for 
T < Tc, (Jf ) oc (1 - T/Tcf^. In other words, in the probe Hmit the critical 
exponent for (Jf) takes the mean-field value 1/2. Does increasing a change the 
critical exponent? Our numerical evidence suggests that the answer is no: for 

1/2 

all a < ac, we appear to find (Jf ) cx (1 — T/Tc) (see Figure 2). 
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Figure 2: The order parameter (Ji), multiplied by K'^/{2n^a^T^), versus the rescaled tem- 
perature T/Tc for different a: a = 0.032 < Oc (green dotted), a = 0.316 < Oc (blue solid) 
and a = 0.447 > Oc (red dashed). The black dot-dashed curve is the function a(l — T/Tc)^^^ 
with a = 160. The green dotted curve is scaled up by a factor of 8 while the red dashed curve 
is scaled down by a factor of 5 such that a, which depends on a, coincides for the green dotted 
and blue solid curves. If we decrease T toward Tc, entering the figure from the right, we see 
that the blue solid and the green dotted curves rise continuously and monotonically from zero 
at T = Tc, signaling a second-order phase transition. The close agreement with the black 
dot-dashed curve suggests that these grow from zero as (1 - TITcfl^. In the a = 0.447 case, 
the red dashed curve becomes multi- valued at T = 1.061 Tc- In this case, at T = Tc, the value 
of Kg(Jj")/(27r'^a^Tc ) jumps from zero to the upper part of the red dashed curve, signaling a 
first-order transition. 




Figure 3: (a) (vr^VT*) versus the rescaled temperature T/Tc for a = 0.316. The 

blue solid curve comes from solutions with (Jf ) = while the red dashed curve comes from 
solutions with nonzero (Jf). For T > Tc, we have only the blue curve, but when T < Tc 
the red dashed curve appears and has the lower (tt'^VTc), indicating a phase transition 

at T = Tc. {-K^VT^) is continuous and differentiable at T = Tc. (b) k|5/ 

versus T/Tc for a = 0.316. The blue solid and red dashed curves have the same meaning as in 
(a). kI^/ (27r*yT^) is continuous but not differentiable at T = Tc, indicating a second-order 
transition. 
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Figure 4: (a) KjH/ (7r*yT*) versus the rescaled temperature T/Tc for a = 0.447. The blue 
solid and red dashed curves have the same meanings as in Figure 3. For T > Tc we have 
only the blue solid curve. At T = 1.061 Tc, the red dashed curve appears and K^n/ (vr'^yT^) 
becomes multi-valued. When T < Tc the red dashed curve has the lowest (n'^VT^), 
indicating a phase transition at T = Tc. Kgf!/ (tt^V^Tc) is continuous but not difFerentiable 
at T = Tc, signaling a first-order transition, (b) k§5/ (2tt'^VT^) versus T/Tc for a = 0.447. 
KgS/ (27r*yTc ) is not continuous at T = Tc, but rather jumps from the blue solid curve to 
the lowest branch of the red dashed curve, indicating a first-order transition. 



As a increases past Uc = 0.365 ± 0.001, we see a qualitative change in 
the thermodynamics. Consider for example a = 0.447. Here again we only 
find solutions with (Jf ) = down to some temperature where two new sets 
of solutions appear, both with nonzero (Jf). In other words, three states are 
available to the system: one with ( Jf ) — and two with nonzero ( Jf ). Figure 2 
shows that as we cool the system, ( Jf ) becomes multi- valued at T = 1.061 Tc- To 
determine which state is thermodynamically preferred, we compute the grand 
potential n. Figure 4 (a) shows nlft/ (tt'^VT^) versus T/Tc. The blue sohd 
curve and the red dashed curve have the same meanings as in Figure 3. We 
immediately see the characteristic "swallowtail" shape of a first-order phase 
transition. If we decrease T, entering the figure along the blue solid curve 
from the right, we reach the temperature T = 1.061 Tc where the new solutions 
appear (as the red dashed curve). The blue solid curve still has the lowest 
K^ri/ [n'^VTc) until T = Tc (by definition). If we continue reducing T below 
Tc, then the red curve has the lowest Kgfi/ (Tr^yT^). The transition is clearly 
first order: K^il/ (tt'^VTc) has a kink at T = Tc. We can also see from the 
entropy that the transition is first order. Figure 4 (b) shows k^S/ (2n'^VTc) 
versus T/Tc. The entropy, like the grand potential, is multi-valued, and jumps 
discontinuously from the blue solid curve to the lowest part of the red dashed 
curve at T = Tc, indicating a first-order transition. 

Notice that a crucial difference between a < ac (second order) and a > ac 
(first order) is that for a > ac the critical temperature Tc is not simply the 
temperature at which (Jf) becomes nonzero. We need more information to 
determine Tc when a > ac, for example we can study fl. 

A good question is: how does increasing a change Tc? Table 1 shows several 
values of a and the associated Tc in units of fixed ^. In the probe limit, a = 0, 
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we have the analytic result from ref. [19] that Tc/ ^l = l/An « 7.96 x 10 ^. For 
finite a, the general trend is that Tc decreases as we increase a. 



a 


Tc/fi X 10^ 





7.96 


0.032 


7.92 


0.316 


4.57 


0.364 


3.62 


0.447 


2.18 



Table 1: Values of Tc/ix, scaled up by a factor of 10^, for various values of a. Recall that the 
critical value of a is Oc = 0.365 ± 0.001. 



5. Discussion eind Outlook 

We studied asymptotically AdS charged black holes in (4 + l)-dimensional 
SU{2) Einstein- Yang-Mills theory with finite a = K^/g, that is, with back- 
reaction of the gauge fields. Our numerical solutions show that, for a given 
value of a, as the temperature decreases the black holes grow vector hair. Via 
AdS/CFT, this process appears as a phase transition to a p-wave superfluid 
state in a strongly-coupled CFT. We have shown that the order of the phase 
transition depends on the value of a: for values below ac = 0.365 ± 0.001, the 
transition is second order, while for larger values the transition is first order. 

As we mentioned in the introduction, intuitively we may think of increasing a 
as increasing the ratio of charged degrees of freedom to total degrees of freedom 
in the CFT. To make that intuition precise, we can consider a specific system. 
One string theory realization of 5/7(2) gauge fields in AdS space is type IIB 
supergravity in (4 -|- 1) -dimensional AdS space (times a five-sphere) plus two 
coincident D7-branes that provide the SU{2) gauge fields [20, 21, 22, 23]. The 
dual field theory is Af = 4 supersymmetric SU (Nc) Yang- Mills theory, in the 
limits of large Nc and large 't Hooft coupling, coupled to a number Nf = 2 
of massless Af = 2 supersymmetric hypermultiplets in the Nc representation of 
SU{Nc), i.e. flavor fields. The global SU{Nf) = SU{2) is an isospin symmetry. 
Translating from gravity to field theory quantities, we have I/k^ oc iV^ and 
l/g"^ oc NfNc, hence a oc y/Nf/Nc, which supports our intuition. We must 
be cautious, however. In the field theory, the probe limit consists of neglecting 
quantum effects due to the flavor fields because these are suppressed by powers 
of Nf/Nc. If Nf/Nc becomes finite, then, for example, in the field theory the 
coupling would run, the dual statement being that in type IIB supergravity the 
dilaton would run, which is an effect absent in our model. We should not draw 
too close an analogy between our simple model and this particular string theory 
system. 

Returning to our simple system, with fully back-reacted solutions we can 
potentially answer many questions: 
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What happens as wc increase a further? As a increases, Tc appears to de- 
crease. An obvious question is whether Tc ever becomes zero. The second-order 
transition occurs because of an instability of the gauge field in the Reissner- 
Nordstrom background [18]. At T = 0, that instabihty only exists for values of 
a below some upper bound [25] . If the transition was always second-order, then 
we would conclude that the transition is only possible for a below the bound: 
if a is greater than the bound, then the instability never appears, even if we 
cool the system to T = 0. If the transition becomes first-order, however, then 
we must rethink the bound: now Reissner-Nordstrom becomes metastable, and 
a phase transition occurs, at temperatures above those where the instability 
appears, so now Tc may go to zero for some a above the bound (if at all). 

Does a new scaling symmetry emerge at zero temperature? The zero- 
tcmperaturc analysis of ref. [25] suggests that the metric exhibits a new, emer- 
gent, scaling symmetry. Our preliminary numerical results suggest that indeed, 
in the approach to zero temperature, an emergent scaling symmetry appears, 
of the form suggested in ref. [25] . The T — )• limit is numerically challenging, 
however, so a firm answer must wait. 

What about the transport properties of the dual conformal fluid, for example 
the electrical conductivity, which at zero temperature should exhibit a "hard 
gap," as explained in ref. [25]? What is the fluid's response to nonzero superfluid 
velocities? In similar systems, sufliciently large superfluid velocities also changed 
the transition from second to first order [31, 32]. What is the speed of sound, 
which need not be the same in all directions since rotational symmetry is broken, 
or the speeds of second and fourth sounds [19, 33, 34]? 

We plan to investigate these and related questions in the future. 
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